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We consider the effect of dynamical localization on the life- 
times of the resonances in open wave-chaotic dielectric cav- 
ities. We show that dynamical localization leads to a log- 
normal distribution of the resonance lifetimes which scales 
with the localization length in excellent agreement with the 
results of numerical calculations for open rough microcavities. 

PACS numbers: 42.55.Sa, 05.45. Mt, 42.25.-p 

The study of lifetime distributions of finite quantum 
systems weakly coupled to a continuum is a subject of 
active experimental and theoretical investigation. The 
nature of the spectrum of resonances depends strongly on 
the nature of the states of the finite system "in isolation" . 
For example if those states are ergodically extended and 
structureless over the system then the resonances will 
show the behavior expected from random matrix the- 
ory, the famous Porter-Thomas distribution in the case 
of a single channel [Q. A close relative of this resonance 
distribution has been measured in quantum dots in the 
Coulomb blockade regime ||^,^ . More recently it has been 
pointed out that optical cavities with partially or fully 
chaotic ray dynamics would have interesting resonance 
properties and efforts have been made to characterize 
their distribution in various limits In a geometry 

which is approximately translationally invariant in one 
direction the wave equation becomes a scalar equation 
with a close formal analogy to the Schrodinger equation 
and the physics of the resonance spectrum becomes essen- 
tially the same for the optical and quantum systems. We 
will henceforth consider cylindrical dielectric resonators 
which are translationally invariant along their axis, but 
can be deformed in their cross-section. The analogue of 
the classical limit of the Schrodinger equation is the limit 
of ray optics when the wavelength of the electromagnetic 
field is much shorter than the typical radius of the cav- 
ity X <ti Rq. We will regularly use the term "quantum" 
to describe properties of the wave solutions which differ 
from the behavior of rays in the same geometry. The mo- 
tion of a light ray within the cavity is identical to that of a 
point mass in a classical billiard and the resulting bound 
states are the analogue of the eigenstates of "quantum 
billiards" Q. However, unless the index of refraction, n, 
is taken infinite, none of these states are truly bound, 
there always being some non-zero probability of escape 
from the cavity. Moreover in the case of a simple dielec- 
tric cavity the escape probability is strongly dependent 



on the angle of incidence of the ray. In particular, rays 
bouncing at the cavity's boundary with an angle of inci- 
dence X smaller than the critical angle, Xc = sin~^(l/n) 
(angles of incidence are defined from the normal to the 
boundary), are transmitted by refraction with high prob- 
ability, while those with x > Xc are trapped by total 
internal reflection, and can only escape with low proba- 
bility by tunneling (evanescent leakage). Semiclassically 
the (dimensionless) angular momentum of the ray in a 
circular cavity is m = nkRa sinx, where fc = 27r/A is the 
wavevector (in vacuum) and Rq is the radius of the cavity. 
Hence a ray with angular momentum m > kRo will be 
strongly trapped whereas one with m < kRo will rapidly 
escape. Correspondingly, resonant states with mean val- 
ues (to) > kRo will have long lifetimes, whereas those 
with mean values less than kRo will have short lifetimes, 
i.e. there is a threshold value rric — kRo for strong escape 
in angular momentum space. In an undeformed (circu- 
lar) cavity to is an integral of motion and there are many 
exponentially long-lived "whispering gallery" resonances 
with TO > kRo. 

For a generically deformed cavity angular momentum 
is not conserved, nor is there any other second constant 
of motion beyond the energy Hence the angular mo- 
mentum can fluctuate. The scale of those fluctuations 
depends on the existence of KAM tori in phase space 
which limit the diffusion in angle of incidence. Beyond 
some critical value of the deformation these barriers are 
destroyed and classical rays with initial angular momenta 
TO much larger than TOc can now diffuse to arbitrarily 
low angular momentum and escape by refraction Q . As 
a result, even for kRo ^ 1 the lifetime of rays start- 
ing with TO > TOc is not exponentially long; the corre- 
sponding level width Tm can be estimated from the dis- 
tance to the critical value in angular momentum space: 
r,„ D/{m — mcY' . (Here D is the effective diffusion 
coefficient in phase space, which in principle can depend 
on TO.) One might then guess that a cavity with such 
chaotic ray dynamics will no longer support any high-Q 
resonances. However this is not necessarily the case, due 
to the phenomenon of "dynamical localization" ||^. It 
is now well-known that just as a random system exhibits 
exponential localization in real-space due to Anderson lo- 
calization, the same kind of destructive interference can 
occur in a chaotic dynamical system, and suppress diffu- 
sion in the relevant phase space pO| |. 

The condition for the onset of dynamical localization is 
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that the diffusion time across the system be longer than 
the Heisenberg time defined by the inverse level spac- 
ing of the cavity: tn ^ hA~^. For longer times than 
tH a wavepacket starts to "resolve" the discreteness of 
the spectrum and the spreading in angular momentum is 
suppressed. Based on an analogy with the kicked rotator 
| pl| , the localization length ^ is determined by the clas- 
sical diffusion rate £), ^ ~ D. Consider a state centered 
around an angular momentum mo such that mo — rric < ^. 
In this case wavepackets can escape before their diffusion 
ceases and the classical picture is adequate. Two differ- 
ent statistical behaviors are possible in this regime. If 
the escape is determined by slow diffusion to a bound- 
ary where escape can occur the survival probability has 
been studied recently by several authors and they 

find characteristic power-law distributions. Here the dif- 
fusion constant satisfies L <^ D L^, and it takes many 
collisions to cross the available phase space (for the op- 
tical cavities the role of the system size L is played by 
2nkRo, the number of angular momentum states avail- 
able). When the diffusion constant is larger, D ^ L^, 
the motion is ballistic in the sense that the phase space 
is crossed in a few collisions; this situation leads to the 
Porter-Thomas distribution of resonance widths and the 
related distributions mentioned above [p],p| jr^ jl^ . How- 
ever when dynamical localization dominates then the life- 
time of a localized state centered around angular mo- 
mentum mo such that mo — mc S> ^ becomes exponen- 
tially longer than the corresponding classical diffusion 
time to the classical emission threshold. Thus one has 
the possibility of high-Q resonances of completely non- 
classical, pseudo-random character, something not con- 
sidered in the optics literature to our knowledge (except 
in a very recent experiment in the microwave regime [ p7| ) . 
It therefore becomes of interest to understand the statis- 
tical distribution of resonance lifetimes in such a situa- 
tion. 

In the localized regime ^/L <^ I, the angular momen- 
tum components of wavefunctions decay exponentially 
away from their centers and one naturally expects ex- 
ponentially long average lifetimes for states centered far 
above the classical emission threshold mo — nic > Re- 
cently Nockel and Stone Q compared the exact lifetimes 
of resonances of quadrupole-deformed microcavities with 
the mean classical diffusion time and found the lifetimes 
to be significantly longer in certain cases; they conjec- 
tured that these discrepancies arose from incipient dy- 
namical localization. Indeed, dynamical localization has 
been shown to occur in certain closed cavities [ p^JTT| 
and a very recent experimental paper confirmed this phe- 
nomenon in microwave cavities of similar shape to those 
studied below jl^. However no detailed study has been 
made of the statistical and scaling properties of the life- 
time in this regime. These are the main topics of the 
current work. Below we will show that this localized 
regime is characterized by a very broad (log-normal) life- 



time distribution with scaling properties directly related 
to the system's localization length 5. 

First, to illustrate the effect of dynamical localization 
on the physical properties of the modes in the open cav- 
ity, we show in Fig.^ the real-space structure of two 
modes of a deformed cylindrical microcavity defined ac- 
cording to the model described immediately below. The 
two modes correspond to exactly the same shape of 
the cavity, corresponding to fully ergodic classical ray 
dynamics, have the same average angular momentum 
^j^2^i/2 0.5nfci?o, but differ in their wavevector and 
as a consequence (see below) differ in their localization 
lengths. As a result one resonance is in the quantum dif- 
fusive regime and the other in the localized regime. The 
qualitative difference is immediately apparent; the dif- 
fusive mode emits much more strongly and has a more 
dense spatial structure due to the large angular momen- 
tum spread in the state. The localized mode on the other 
hand emits weakly and appears to have a caustic similar 
to a regular whispering gallery mode, but a closer look 
at its spatial structure shows the pattern of nodes has 
an irregular character entirely different from the usual 
whispering gallery modes of circular resonators as can be 
seen on Fig. ^. 

We now define the model corresponding to Figure 
1 and 2. We consider an optically inactive, cylindri- 
cal microcavity with an index of refraction n > 1. 
The cross-section perpendicular to the cylinder's axis 
is given by a circle perturbed by M harmonics of ran- 
dom amplitude < a£ < l/l, R{(t>) = ^o[l + 
^^f^2 '^^ cos(£0)]. The average roughness of the surface 
is defined as k = < n'^^cj)) >0, k(0) = {dR/ dcj)) / Rq. 
This model was introduced by Frahm and Shepelyan- 
sky with the condition of perfectly reflecting walls, 
and they referred to it as the rough billiard to contrast 
with the smoother quadrupolar deformations considered 
by Nockel and Stone Q . However the spatial wavelength 
of the roughness is still assumed to be large compared 
to the wavelength of the resonance. The advantage of a 
rough boundary is that the transition to classical chaos 
is achieved with much smaller amplitude of deformation 
making it easier to explore the parameter regime of fully 
chaotic classical motion and dynamically localized " quan- 
tum" behavior. As we shall see below, the open rough 
billiard has scaling and statistical properties essentially 
identical to a quasi-one-dimensional disordered system, 
whereas the quadrupole billiard does not. For the rough 
billiard the classical dynamics can be well approximated 
by a discrete map for which Chirikov's overlap criterion 
|l9| gives an estimate of the critical roughness Rc above 
which the classical dynamics becomes fully chaotic as 
Kc ^ M~^/'^. The two deformation parameters M and 
K allow to reach a classically fully ergodic regime char- 
acterized by a diffusion constant (averaged over angular 
momentum) D = ^{RnkRo)^ for R ^ Rc and quantum 
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FIG. 1. Top : Intensity plot of a resonance with nkRo = 50 
in the rough cavity with n = 0.08, M = 15 and n = 2.5. Red 
color corresponds to the maximum of the intensity, and blue 
to the minimum. Bottom : Intensity plot of a resonance with 
nkRo ~ 150 for the same set of parameters as above. 

mechanically, one gets a localization length ^ ~ 13 ^ so 
that the dynamically localized regime is determined by 
K^nkRo <^ 1 [0. Keeping parameters of the cavity fixed 
and varying kRo one is able to access states with greatly 
diff'erent localization lengths. 

We restrict ourselves to the simplest case of TM- 
polarized electric field E{r) parallel to the cylinder axis 
for which both the field and its derivative are continuous 
at the cavity's boundary. This restriction is primarily 
for convenience, the TE modes obey a slightly different 
scalar equation which can be treated in a similar man- 
ner. It should be mentioned however that semiconductor 
quantum cascade micro-cylinder lasers studied in |^ emit 
solely in the TM mode due to a selection rule. Maxwell's 
equation reduce then to a single scalar wave equation, 



FIG. 2. Details of the wave intensity corresponding to the 
top resonance shown in Fig. ^ The intricate structure of the 
wave intensity is due to classically chaotic boundary scatter- 
ing and makes the resonance clearly different from a standard 
whispering-gallery like high-Q resonance. 

c~^dfE{r,t) = n^(r, 0)Vr£'(r, i), where the refraction 
index satisfies n(r) — n inside the cavity, and n(r) — 1 
outside. 

We use the approach in which the resonances widths 
in wavevector are given by the imaginary part of the 
wavevector of the quasibound states defined by the fol- 
lowing matching conditions. First we expand the electric 
field in the angular momentum basis (r = |r|) 

E{v,t)^e-''^' ^ i^vl™ (fcr) e*™*, (1) 

m— — oo 

where 

. _ J a,nH+ (nkr) + /3„ii?„ {nkr) if r < R{(j>), , , 

" ~ 1 l,nH+ (kr) otherwise. ^ ' 

This expansion corresponds to the so-called Siegert 
boundary conditions po|] in which the states have only 
an outgoing component at infinity {H^{x) are Hankel 
functions of first and second type, respectively). Such 
boundary conditions cannot be satisfied for real k and 
are only satisfied for discrete complex k. It can be shown 
that the imaginary part of the k values which satisfy Eq. 
(2) are the poles of the true unitary (on-shell) S-matrix of 
the scattering problem. From the expansion coefficients 
in Eq. (^) we define vectors a, /3 and 7. The fields inside 
and outside the cavity are related by the continuity of the 
field and its derivative on the boundary and (after inte- 
gration around the boundary) one of these equations can 
be used to eliminate 7 leaving a linear relation between 
a and /3. The matrix expressing this relation we call S. 
Moreover the regularity of the field at the origin r = 
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implies a = /?, thus a secular equation for the resonant 
values of k is obtained of the form: 



(3) 



We use the notation iS-matrix because in the case of 
a closed billiard the matrix so-defined is actually the 
unitary S-matrix of the scattering problem of a wave 
incident outside the impenetrable billiard |^.^. In 
our case the matrix so-defined is non-unitary for any 
n < oo and for real k, the eigenvalues of S have the form 
Xr = exp(i((/3r + iSr)), where both (pr and > are 
real functions of momentum k. The subscript r numbers 
states in a deformed cavity where angular momentum is 
not conserved. Exact quantization of the cavity - solving 
Eq. (^) exactly - implies = 1 so that the exact imple- 
mentation of this procedure requires finding a complex 
k — q — i"f such that </?r(<? — *7) = ^ril — i"f) — 0. The 
corresponding inverse lifetime Tr is then given by c times 
7 via the dispersion relation, Tr = cy. However approxi- 
mate lifetimes can be found by a much more efficient pro- 
cedure which is extremely helpful if one wishes to study 
full distributions as we do here. First, it is straightfor- 



ward to show 1 23 1 that when 7 is small, simply finding 
the complex eigenphases = exp(i((^r + *<5r)) for real k 
determines the imaginary part of k on resonance by the 
relation: 7 = —S/(p' where ip' is the derivative of the real 
part of the phase with respect to momentum for real k. 
Moreover, since this derivative can be shown to be slowly- 
varying on the scale of the level spacing Afc, it is not 
necessary even to quantize the real part of the phase (i.e. 
to find the real k which makes ip{k) = 27r x integer). The 
function ip' can be easily calculated for the circular cylin- 
der and this relation and the assumption of slow variation 
of the derivative can be confirmed explicitly (for the case 
where 7 is small). Therefore we can generate large life- 
time ensembles simply by diagonalizing the matrix S for 
real k and extracting the imaginary phase d, by which 
means we generate ~ 2nkRQ inverse lifetimes per diago- 
nalization. This procedure is motivated by the work of 
Doron and Frischat who noted that the statistical proper- 
ties of closed billiards changed little away from the exact 
quantization condition (in their case it was the distri- 
bution of splittings of semiclassically degenerate states) 
. The linear relation between Sr and 7^ has been in- 
dependently proposed earlier by Hackenbroich ]25t and 
demonstrated for the case of the circle. 

To best relate the localization properties of the eigen- 
states, which apply to a closed cavity, to the distribution 
of lifetimes in an open cavity, we employ a perturbative 
formalism which was recently developed specifically to 
treat open optical resonators ||^ (it is similar in spirit to 
well-known quantum perturbative scattering approaches 
such as R-matrix theory in the single-level approxima- 
tion). According to that theory narrow resonance widths 
r <C A (A is the resonance spacing) can be computed 



from the expectation value of an antihermitian operator 
V taken over eigenstates | a^"^) of the matrix Ai, which 
describes some effective "closed cavity" 

r = 6o(a(°) 11^1 «(")) =boJ2 (4) 



Explicitely, 

M = [J' J') - {l/n){j'H+') {H-H+y' {H-J) , (5) 

V is the antihermitian part of M, and the matrix ele- 
ments (ZZ) are defined as 



2tt 



Zi (fc) Zm (k) e 



(6) 



Z,n{k) and Zmik) stand for either H^{kR{(j))), the 
Bessel function Jm {nkR{(j))), or their derivative. The 
coefficient bo in (^ depends only on the hermitian 
part Ho of ^A, i.e. it is determined by the properties 
of the "closed system", bo^ = (a(°) |aHo/5fc| 



fO)* 



{J'J")r 



„(0) 



and can be regarded as a nor- 



malization factor. An eigenstate, localized at angular 
momentum mp S> mc ~ kRo, will have exponentially 
small width F. Therefore, for such a resonance in the 
semiclassical limit F <C A, and Eq. (Q) is appropriate. 
More details on this formalism can be found in Q . 
For an exponentially localized state one generally has 
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I'm - mo\ 

2e 



(7) 



We consider the regime of large localization length 
^ ~ R^k^Ro ^ 1. Since the rough billiard is classically 
chaotic, the phases of the coefficients change rapidly 
with the angular momentum, and it is natural to assume 
that its correlation function satisfies 



(8) 



where the average is performed over an angular mo- 
mentum interval I € [mo — S£/2, mo + S£/2] such that 
1 Si < This behavior is illustrated in the inset to 
Fig. ^ for one typical set of cavity parameters, corrobo- 
rating the validity of the assumption (||) . 

As follows from Eq. and the definition of V, the 
matrix elements Vmm' in angular momentum representa- 
tion vary on a scale of Sm ~ RkRo S> 1 , and therefore we 
can replace the product aj^am' in Eq. (^ by its average 
value (Q;*jQ!m') over the interval |to — m'l ~ RkRo- To- 
gether with (p|) this leads to the diagonal approximation 



nkRo 

{ao \V\ ao) « ^ I 

\m\—kRQ 



(9) 



The matrix element 
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c.c. 



(10) 



includes both the refractive (classical) escape from the 
resonator (for \£\, \£'\ < rric = kRo) and the "tunnel- 
ing escape " (corresponding to evanescent leakage , for 
l-^l, l^'l > trie). To evaluate the sum over angular mo- 
menta we use the stationary phase-based technique, de- 
veloped in in the context of the calculation of level 
splittings in rough billiard. The "classical" refraction con- 
tribution is found to be 



mo >mc 



k^Rp 



■ exp 



mp — kRo 
I 



(11) 



This result shows that an exponentially-long lifetime can 
be due to the exponentially-small wavefunction com- 
ponent leaking outside the classically totally-internally- 
reflected region. To see if this process controls the life- 
time we need to compare this result with the direct "tun- 
neling escape" contribution to the lifetime. The latter 
process involves angular momenta only above emission 
threshold rric. For an estimate, it is then sufficient to 
evaluate the linewidth of the state above rric in the circu- 
lar cavity, which can be thought of as a state with zero 
localization length. We find that the tunneling contribu- 
tion is also exponentially small (ti 3> 1), 

Co">m. = -(l/n)Im [H+^^_,ikRo)/H+^{kRo)] 

~ exp (^2^ tuq — ml + 2too In(mc) 

-2mo ln(mo + \J ml ~ m^)^ . (12) 

Competition between the classical escape [Eq.(|ll|)] and 
the tunneling [Eq.(12)] is strongest for (mo — mc)/mc <C 
1, i.e. when the width of the tunneling barrier is smallest. 
Comparing the two contributions in this region we find 
that r'^i^'^'^ > r*""" for C > VkTk,, and T*""" > T'^''^"" in 
the opposite limit. Restriction on the range of ^ weak- 
ens as one moves to higher angular momentum, and for 
(too — rric) I rric ~ 1 the classical escape mechanism al- 
ways dominates over the tunneling one. The tunneling es- 
cape therefore is relevant only for very small deformations 
K <JC {kRo) , which produce short localization length. 
Thus, lifetime of the states with very short (essentially 
zero) localization length is determined by the tunneling 
escape, whereas that of the more extended states with 
^ ^ 1 by the classical emission from their exponentially 
weak tails at rric — kRo- 

Having established a relation between the lifetime and 
the localization length of a corresponding closed cavity, 
the distribution of lifetimes then follows from the distri- 
bution of localization lengths. In one-dimensional and 
quasi-one-dimensional disordered systems it is now well- 
established that the distribution of the inverse localiza- 
tion length is typically normally distributed around an 
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FIG. 3. Distribution of the inverse participation ratio p for 
n = 3, M = 15, and {nkRo,R) = (150,0.02) (solid line), 
(100,0.03) (dashed line) and (50,0.06) (long-dashed line). 
RnkR has been kept constant, which result in a stable av- 
erage IPR (p) = 0.1 ± 0.01. Inset: Normalized correlation 
function (^) for nkRo = 100, and k — 0.03. 

average l/^o (2^]- Moreover Frahm and Shepelyansky 
have explicitly shown pT| that the problem of the rough 
billiard maps onto a variant of the kicked rotor problem 
and hence to an ensemble of band random matrices , 
which also describe quasi-ld disordered systems. Here 
the angular momentum index plays the role of the site 
coordinates in the disordered systems, with an ideal lead 
(the continuum) accessible at rric = kRo. Collisions with 
the rough boundary correspond to random hopping be- 
tween sites, which are at most RnkRo lattice spacings 
apart [ |TT[ . Once the critical angular momentum, corre- 
sponding to the classical emission border rric, is reached, 
the wavepacket escapes the system and never returns. 
Hence we may assume that the inverse localization length 
in our problem has an approximately normal distribution 
around its mean (our ensemble here is of boundary real- 
izations) 



P{0 exp[- 



(l/e-l/Co)^ 
2a2 



(13) 



Therefore, as follows from (|ll|) and (|13|), the resonance 
lifetime is distributed log-normally, 

PmHt)) - exp[- J/^'^^^X^^ ,, ]d{Ht)) (14) 
2(7^ (too — kRo)^ 



where 



In(to) = (too - kRo)/ 5,0 



(15) 



and the derivation is done in a leading logarithm approx- 
imation, so that the pre-exponential factor in ( pT| ) is ne- 
glected. This result is then very natural: the distribution 
of lifetimes in open dynamically-localized cavities is log- 
normal for the resonances localized far from the classical 
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emission threshold mo ^ fci?o- This is entirely analo- 
gous to the conductance distribution of localized chains, 
which will be log-normal for a fixed distance from the 
ends (see for log-normal distribution of delay 
times/resonance widths). We note that the relationship 
between dynamical localization and Anderson localiza- 
tion was first placed on firm footing in a seminal paper 
by Fishman, Grempel and Prange px| ]. 

Equation essentially relies on the two assump- 
tions: first , Eq. (||) that the phase of the wavefunc- 
tion components are randomly distributed with no long- 
range correlations and second, that the eigenstates are 
exponentially localized with a normal distribution of lo- 
calization lengths. We now test the validity of these two 
assumptions for our rough microcavities. The validity of 
Eq. (^) is confirmed by the sharp drop of the correlation 
function {al_^^ai)t/ {\at\ )i for m > which is clearly 
seen from the numerical results presented in the inset to 
Fig. ^. The localization properties are investigated by 
computing the distribution of the Inverse Participation 
Ratio (IPR) defined as p = J2m l^mW Em la^np. The 
IPR measures the inverse number of effective eigenstates 
components and thus allows to distinguish between lo- 
calized and delocalized states. Generally, in the local- 
ized phase (p) is independent of the system size since at 
most <^ L sites contribute to the sum and (p) ~ 
In the other limit of ergodic states, all sites contribute 
equally and (p) ~ in this case. Between these two 
limits, a variety of behaviors may occur depending on 
the inner structure of the eigenstates. In Fig. || we 
show IPR distributions for three different parameters sets 
corresponding to the same average localization length 
^ ^ [RnkRo)^ = 9. The three distributions are indeed 
stable under parameter variations keeping ^ constant and 
this shows that not only the average IPR/localization 
length but also the full IPR distribution obeys a 
one-parameter scaling with RnkR. The situations is very 
similar to the one studied in Ref. ||2^ for BRM with a 
band width 1 <C 6 ~ ^f2£, for which the IPR distribution 
was analytically computed. Similar deviations from ||29| ] 
as seen on Fig. ^ are also present for BRM with not too 
large band widths, so that these numerical results con- 
firm the universality of the dynamically localized regime, 
quite analogous to quasi-one-dimensional disordered sys- 
tems. We also illustrate this exponential localization by 
showing one typical state in the inset to Fig. |[ 

Having tested the validity of the main assumptions 
on which ( p^ ) relies, we present in Fig. ^ distribution 
of lifetimes for the classically-chaotic, dynamically local- 
ized regime of the open rough microcavity. The distribu- 
tions shown correspond to resonances centered in inter- 
vals of width 5m/ (nkRo) — 0.1 around angular momenta 
TOo/(nfci?o) = 0.5, 0.6, 0.7 and 0.8, well above the clas- 
sical threshold mc/ (nkRo) « 0.29. Clearly, the distri- 
butions are log-normal and their widths increase as one 
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FIG. 4. Lifetime distribution for n = 3.5, nkR « 110, 
K = 0.018 and Af — 15 for states with average angular mo- 
mentum mo/{nkRo) — 0.5 (circles), 0.6 (squares), 0.7 (dia- 
monds) and 0.8 (triangles). This last distribution is fitted 
according to Eq. (|lj) with a ^ 0.04 (dashed line). The 
emission border is at rric ~ 0.29nkR. Each distribution is 
constructed from 8000 to 13000 lifetimes obtained from 1013 
boundary realizations of the rough cavity. Inset : Typical lo- 
calized eigenstate for the same parameter set. The dashed line 
indicates an exponential decay corresponding to a localization 
length of 5 = 16. 



moves away from the classical emission border nic = kR^. 
Furthermore, the agreement with Eq. (^4|) is quantita- 
tively confirmed by a direct fit of the broadest of these 
distributions (see dashed line on Fig.^). 

We expect by analogy to the scaling theory of local- 
ization that the logarithmic average of the lifetime will 
exhibit a universal scaling behavior. This expectation is 
confirmed by the data shown in Fig. |^ where we present 
numerical results for the scaling obeyed by Into- Log- 
averaged lifetimes for different parameter sets have been 
computed for at least 2000 lifetimes in narrow energy 
windows dm/ (nkRo) — 0.2 around given angular mo- 
mentum mo > mc for different values of indices of re- 
fraction 1.5 < n < 4, wavelengths 75 < nkRo < 180 
and roughnesses k. All presented results are in the local- 
ized regimes ^ < mo — mc and the corresponding curves 
have been put on top of each other by a one-parameter 
scaling. Fig. ^ demonstrates the validity of the linear 
relation (^5|) as is indicated by the straight line. The 
exact parameter dependence of the scaling can be de- 
duced from the analogy to the scaling theory of localiza- 
tion. In this case the Thouless conductance g = F/Ae 
is the scaling quantity (or its logarithm in the localized 
regime); here F is the resonance width, and Ae is the 
mean level spacing. In our case T = cy (where 7 is 
the width in momentum space), but Ae differs from the 
corresponding Schrodinger equation, due to the differ- 
ent dispersion relation for the wave equation. Taking 
this into account one finds that the analogue of the di- 
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mensionless conductance is g ^ n^kcjRQ. and it is the 
logarithm of this dimensionless quantity which we plot 
against (mg — kR)/S^* . Fig. ^ allows to identify the scal- 
ing parameter ^* with the localization length up to a free 
parameter. That this scaling holds for the rough cavity 
demonstrates that localization length is independent on 
the angular momentum as is expected for an homoge- 
neously diffusive system. The situation is fundamentally 
different for a quadrupolar cavity (M — 2) as can be seen 
on Fig.|^ (see black diamonds). Obviously one scaling pa- 
rameter is not sufficient to bring the corresponding curve 
on top of the other ones satisfying ( p^ . This indicates an 
angular momentum dependent diffusion constant which 
directly follows from the effective local map derived for 
this particular case |3|] . Furthermore, in the regime cor- 
responding to the presented data for the quadrupolar de- 
formation, small invariant torii and islands of stability 
still survive, resulting in strongly localized wavefunctions 
with short localization length determined essentially by 
the size of remaining classical structures. Because of this, 
and unlike the situation in the rough cavity, lifetime of 
such states is determined by the tunneling escape (see 
discussion after Eq.(|l2|)). Therefore a clean demonstra- 
tion of exponential dynamical localization is difficult in 
the quadrupolar billiard. 

Further confirmation that the extracted scaling param- 
eter is indeed related to the system's localization prop- 
erties is given in the inset to Fig. || where ^* is plotted 
against the diffusion constant D — [RnkR^)^ as derived 
from the effective rough map ||ri| . This inset gives an un- 
ambiguous confirmation of the above derived relation be- 
tween localization length and log-averaged lifetime. Note 
that ^* has a linear dependence on the diffusion constant 
D even at small D where the relation £^ ^ D does not 
hold. For large D however, the relation £,* ^ £, holds. 

To summarize, we have presented the first study of 
the lifetime distributions of a quantum-chaotic open sys- 
tem in the dynamically localized regime. This study was 
greatly expedited by the linear relation between the com- 
plex phases of the eigenvalues of the nonunitary scatter- 
ing matrix away from exact quantization and the imag- 
inary part of the corresponding exactly quantized com- 
plex wavevector, which has allowed us to generate suf- 
ficiently large lifetime statistics to demonstrate the log- 
normal form of their distribution. The lifetime distri- 
bution has been derived analytically assuming a normal 
distribution of inverse localization lengths and phase ran- 
domness of the wavefunctions, using a recently developed 
semiclassical method, the usefulness of which is thus fur- 
ther demonstrated and confirmed numerically. The 
log-normal distribution is a hallmark of localized disor- 
dered systems and hence our results deepen the analogy 
between dynamical and Anderson localization and point 
out an optical observable which can in principle be mea- 
sured to demonstrate this distribution. The possibility 
of high-Q resonances in deformed rough cavities (which 
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(m„-kR)/^ 

FIG. 5. Scaling of the average logarithm of lifetime (Inio) 
vs. (mo — kR)/^* for parameters 95 < nkR < 195, 
2.5 < n < 5, 0.008 < R < 0.03 and 10 < M < 30. Open 
symbols correspond to the rough deformation for which a 
one-parameter scaling exists. The full diamonds correspond 
to a quadrupolar billiard (M = 2) for which the diffusion is 
affected by classical invariant structures in phase space. Inset 
: Localization length as extracted from the scaling shown on 
the main pannel. 

are nonetheless smooth on the scale of the wavelength) 
should be of interest in optical studies of scattering from 
small particles, however their random nature seem to 
make such resonances unsuitable for applications. 
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